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ABSTRACT 

Recent observations suggest the incidence of strong intervening Mgn absorption sys- 
tems along the line-of-sight to gamma ray burst (GRB) afterglows is significantly 
higher than expected from analogous quasar sightlines. One possible explanation is a 
geometric effect, arising because Mgn absorbers only partially cover the quasar contin- 
uum regions, in which case Mgn absorbers must be considerably smaller than previous 
estimates. We investigate the production of abnormal absorption profiles by partial 
coverage and conclude that the lack of any known anomalous profiles in observed 
systems, whilst constraining, cannot on its own rule out patchy Mgn absorbers. 

In a separate test, wc look for differences in the distribution function of Mg n equiv- 
alent widths over quasar continuum regions and C in] emission lines. We show that 
these anomalies should be observable in any scenario where Mgn absorbers are very 
small, but they are not present in the data. We conclude that models invoking small 
Mgn cloudlets to explain the excess of absorbers seen towards GRBs are ruled out. 
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1 INTRODUCTION 

The recent claim that strong Mgn absorption along the 
sightlines to gamma ray bursts (GRBs) is sig nificantly more 
likely than along the sightlines to quasars (jProchter et alj 



likely 

l200d . hereafter P06) is puzzling. Subsequent results sug- 
gest that Mgn absorbers are unique in exhibiting this trend, 
C IV absorbers displaying population stat istics independent 
of the nature of the b ackground source (|Teios et al.ll2007l ; 
ISudilovskv et al.ll2007l) . One possibility is that some of the 
Mgn absorbers seen towards GRBs are intrinsic, although 
the high velocity separation between GRB and absorber 
coupled with the low velocity dispersion in the observed 
Mgn components prove a challenge to theoretical models. 
Therefore, it remains of interest to investigate possible ex- 
planations for the effect in terms of intervening absorbers. 

This le t ter w ill focus on the model proposed by 
iFrank et all (|2006l , hereafter F06), which suggests geomet- 
ric effects are responsible for the discrepancy. In this picture, 
the Mgn absorbers are somewhat smaller than the beam size 
of the quasar continuum region, whereas GRB beam sizes 
are smaller still. This would evidently lead to a difference in 
the observed distribution of absorber properties. 

However, such small absorber sizes do not fit well with 
the established observational picture. Information on ab- 
sorber sizes comes from studies of galaxies near quasar ab- 
sorber sightlines, and from comparing absorber and galaxy 
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number densities. These give typical radii ~40kpc for 
Mgn complexes, which may contain several velocit y com- 
ponents, for a rest equiv alent width limit of 0. 3A jsteidell 
1 19931 ; IZibetti et alll2007n . For l ower equivalent width lim- 
its the inferred sizes are larger (|Churchill et al.lll999l ), and 
projected Mgn absorber-galaxy correlation f unctions show 
a sig nificant excess out to over 1 Mpc (e.g. iBouche et al.1 
2004). For individual velocity compone nts of low ionizatio n 
species, sizes are less than ~300-400pc (|Rauch et al.ll20o3 ). 

On the other hand, quasar continuum regions are gen- 
erally considered to be orders of magnitude smaller. Impor- 
tant for the investigation presented here is the significant 
difference in the physical size of the ultraviolet continuum 
emitting region and the broad emission line region in lumi- 
nous quasars. Constraints of the physical size of the broad 
line and co ntinuum regions are pro vided by reverberation 



studies (e.g. Kaspi et al.| | 2005 , 120071 ) and gravitational lens- 
ing studies ( Lewis et alJll99S )■ Recently the latter method 
was used by Wavth et al.l ( 20051 ) . giving values of 0.02 pc for 
the continuum region, with a broad-line region > 3 times 
larger (possibly considerably larger than this lower limit). 

A comprehensive analysis of possible routes to solv- 
ing the GRB discre pancy has recently been presented by 
iPorciani et al.l (|2007| ) , wherein it is suggested that the model 
in F06 has severe weaknesses. These relate to the neglect 
of multiple clouds in the F06 model, and the difficulty in 
producing dou blet ratios which are consistent with obser- 
vations. Whilst IPorciani et all l|2007h propose strong objec- 
tions to the model, they do not rule it out. Furthermore, the 
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idea appears to have been gaining po pularity and has been 
used to interpret recent observations l|Hao et al.ll2007l ). 

In this letter, we provide direct observational evidence 
against the F06 model. First, we briefly discuss whether 
patchy absorbers can be ruled out on the basis of their 
absorption profiles and line ratios alone (Section [2}, and 
conclude that this is not always the case. We then use the 
differing size of the quasar continuum and C ill] emission re- 
gion to provide a definitive test for the F06 model (Section 
[3]) . We provide a summary in Section [4] 



2 PROFILES OF SINGLE VELOCITY 
COMPONENT PARTIAL ABSORBERS 

Suppose quasars emit light with a hard-edged circular beam, 
radius R e . Then a circular, single component Mgn absorber 
can be parametrized by its central column density (A^Mgii), 
velocity width (v) and its radius as a fraction of the beam 
radius (Ro/R e ). In practice the 'beam width' is an angular 
line-of-sight effect and will scale accordingly as Ro = 8Da(z) 
where z is the redshift of the absorber and Da is the angu- 
lar diameter distance. However, this geometric effect will be 
neglected as Da (z a bs) — -DA(zem). We set R e = 1 without 
loss of generality. 

For Ro < 1, providing that spectra of adequate S/N 
and velocity resolution are obtained, saturated profiles with 
a systematic offset above the zero-flux level (resulting from 
light 'leaking' around the edges) will be seen. However, it is 
not realistic to model absorbers with a hard edge; there must 
be a region around their edges where the column density 
drops. In the absence of any favoured theoretical or empirical 
model, we write: 



N- 



Mgll 



(R) = 



No 



R < Ro 
R > Ro 



(1) 



During the drafting of this paper, a revised version of P06 
independently proposed a similar model, using a three di- 
mensional distribution: 



p{r) 



Po 
Po 



r < r 
r > ro 



(2) 



This latter model is less convenient for our purposes of com- 
puting high resolution profiles, something which P06 did not 
attempt. However, one may roughly relate the two distribu- 
tions simply by integrating the second to yield n — k — 1 
outside the core radius. 

The resulting observed profile is calculated numerically 

as: 



F(A)oc [ 2wRe- Tx{NM ^ l( - R) ' v) dR 
Jo 



(3) 



where t\ (N, v) is a standard Voigt optical depth profile for 
Mgn A2796, and F(A) is the observed flux. We separately 
calculate the somewhat weaker profile Mgn A2803 so that 
we can fit for both components simultaneously. We take 
v = 5 kms -1 , which is the me an broadening for individ- 
ual Mgn clouds (|Churchillll 19971 ). 

The model spectra are realized at high resolution (SX = 
0.01 A in the absorber restframe). Gaussian random noise 
is added with S/N=100, to approximate the best available 
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Figure 1. The reduced x 2 of the best fit pure Voigt profile to a 
partial absorber with a core radius Ro (R e = 1), outside of which 
the column density decays as R~ 1 (dashed lines) and R~ 2 (solid 
lines; the shaded regions indicate the variance of \ 2 of individual 
fits, which is the same for the different n values but differs accord- 
ing to line strength because of the different number of pixels used 
in constraining the fit). The red, green and blue lines indicate the 
results for observed column densities of log 10 N ^ s = 12, 14, 16 
respectively. Note this means the intrinsic central column den- 
sity increases with decreasing radius. As an example, the inset 
figure shows the profile for Rr 2 law, iV D bs = 10 14 , Ro/Re = 0.5, 
with the best fit plotted in green. By eye the partial coverage is 
virtually impossible to detect. 



echelle spectra. We then fit a Voigt profile to the two com- 
ponents, assuming the correct ratio of oscillator strengths. 
The resulting reduced \ 2 1 measured over the region where 
the departure from continuum level is significant at the la 
level, gives an indication of how well the partial absorber can 
fool one into believing the feature is due to a full coverage 
absorber, including the effect of the doublet ratio. 

We consider a grid of models of varying JV b s and Ro, 
where JV b s is the observed (fitted) column density, not the 
intrinsic column density of the partial absorber. This is 
achieved by increasing, for each model, the intrinsic column 
density No until the best fit gives the target N b s . The pro- 
cedure ensures each case we consider is a candidate for an 
observable system, and automatically includes cases where 
iVo is considerably larger than is typically seen in Mgn ab- 
sorbers but Ro is very small, which should be the most chal- 
lenging case to disguise. 

The results for n = 1,2 (chosen for the reasons given 
below) are shown in Figure 1 by dashed and solid lines re- 
spectively with log 10 iVobs = 12, 14, 16 in red, green and blue. 
The shaded regions indicate the ler variance on the reduced 
X 2 for individual fits: note that this depends on the number 
of pixels used in its construction, and so will depend on iV bs 
(but not significantly on n or Ro, since iV bs is held constant 
as these are varied). 

F06 reproduce their geometric results for profiles as soft 
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as n = 1 (k — 2) with core radii Rq ~ 0.25. According to our 
investigations, individual absorbers in this scenario would be 
observationally indistinguishable from the complete cover- 
age case, with reduced x 2 ~ 1 for all fits down to Rq — 0.01. 
As one increases n (making the absorber edges harder), it 
rapidly becomes easier to distinguish the cases; for example, 
for n = 2 and 7?o ~ 0.25, the high resolution spectra sim- 
ulated here will highlight a problem. Nevertheless, the \ 2 
values presented in Figure f represent the poorest possible 
quality of fit given that the model spectra are of excep- 
tionally high quality and only allow for fitting one velocity 
component. 

Given the results described above, it may appear possi- 
ble to have a situation in which partial absorption is present 
but has not been detected in the high resolution data. The 
physical explanation is simply that the shape of the line de- 
pends on N only logarithmically, so allowing N to vary as 
a shallow powerlaw across the beam does not have a very 
strong effect: nonetheless, according to F06, one obtains sta- 
tistical variations in the equivalent width distribution ac- 
cording to the size of the background quasar. 

However, with such very slow powerlaw decays, there 
are further problems to be investigated. On the one hand, 
the F06 statistical result is valid only in a regime where the 
differing equivalent widths of observed absorbers is gener- 
ated by the differing radii of intersection through a popu- 
lation of very similar clouds. On the other hand, there is 
a very large area over which low column densities may be 
observed on the periphery of such systems. In particular the 
observed column density distribution of a single system may 
be numerically calculated to scale as dP(W)/dW ~ W~ :i . It 
seems unlikely that one could contrive a cosmological distri- 
bution of absorbers which provide a rescue from such a steep 
powerlaw to give th e observed dP(W)/dW ~ W -1 over- 
all di stribution (e.g. iNaravanan et al" 2007; IChurchuT et al.l 
1 19991 ) whilst preserving the F06 result for n = 1. 

In conclusion, it may be that soft edges can be invoked 
to alleviate somewhat the high resolution fitting problems 
associated with patchy absorption. However, whilst the limit 
of extremely shallow outer slopes yields very convincing 
Voigt profiles, one obtains a new set of statistical challenges. 
For the rest of this letter, we will assume sharper cut-offs 
which do not suffer from these problems. 



3 Mgn ABSORBERS CLOSE TO Cm] AI909 
EMISSION IN QUASAR SPECTRA 

If, as in F06, Mgll systems are somewhat smaller than 
quasar continuum beam widths, it follows that they can- 
not cover the broad emission line region (Section [TJ, except 
with a much reduced column density which must be neg- 
ligible according to the statistical argument at the end of 
Section [2] Accordingly, one can approximately model the 
spectrum over a broad emission line as: 



over the continuum region only) will be seen with equivalent 
width W: 



F(X) = C(A)e" TM s Il(A) + L(A) 



(4) 



where F is the received flux, C is the continuum emission 
flux, L is the broad-line emission flux and TMgll is the optical 
depth due to an intervening Mgll system. Assuming Q, 
an absorber with "intrinsic" equivalent width W (as seen 



W' = 



w 



1 + L/C 



(5) 



where L(X) and C(A) are assumed to be constant over the 
absorption feature (compare typical widths of 20kms~ 1 and 
5000 kms for a Mgll component and quasar broad line 
emission respectively). This result applies for any absorber 
size that is smaller than the continuum region. Possible com- 
plications are discussed in Sections 13. 3l fe lX4l but essentially 
this simple model is adequate for our investigation. Equation 
(JS| implies that, for small absorbers, the observed number 
of absorbers above a specified equivalent width threshold 
may be reduced significantly over broad line emission fea- 
tures because of the 'diluting' effect of the light from the 
large broad emission line region. It is thus possible to apply 
a test similar to the ge ometric intrinsic /extr insic test of the 
Lya forest described in lSargent et al.l l|l980f ). 



3.1 Observational Sample 

The Mgll AA2796,2803 absorber sample of IWild et all 
|2007l . Section 2) is used as the basis for the statistical in- 
vestigation of the observed absorber properties in the vicin- 
ity of the C ill] A1909 emission line. It is derived from a 
S /N-limited search of 32 278 quasar spectra contained in 
the Sloan Digital Sky Survey (SP SS) Data Release 4 (DR4) 
l|Adelman-McCarthv et aLll2006l ). The sample includes ab- 
sorbers down to a restframe equivalent width (W) of 0.5 A 
over an extended redshift range. For our purposes we con- 
fine attention to the 2 831 absorbers with Wq 796 SsO.9 A and 
— 180 < <5A/A < 180 where <5A is the wavelength shift be- 
tween the emission line and absorption line in the quasar 
restframe: 



8X = 



2796A(1 + ZMgii) 



r) 



1909 A 



(6) 



Here ZMgtt is the redshift of the absorber and 2 qU asar is the 
redshift of the quasar. The ±180 A wavelength interval was 
chosen to provide sufficient absorption lines in the compar- 
ison sample. The W^ 2796 ^0.9 A limit was chosen to provide 
an absorber sample consisting of a large number of strong 
Mgll absorbers but our conclusions are insensitive both to 
the precise Wo-limit employed or the wavelength extent used 
to define the comparison sample. 

For the test below, we need to predict the absorbers as a 
function of SX, without the correction described by equation 
(JS| which will be considered separately. The expected num- 
ber will be a function of (a) the accessible redshift path, (b) 
the mean redshift of the absorbers and (c) the S/N in the 
quasar spectra. The accessible redshift path as a function of 
S X is readily calculated using the redshifts of the base sample 
of 32 278 quasars. The small systematic change in the num- 
ber of absorbers due to evolution with redshift can be cal- 
culated from the mean absorber redshift as a function of 5X, 
combined with the absorber number density as a function 
of red shift, oc (1 + ZMgii) ' 23 , determined by iNestor et al.l 
(2005). The probability of detecting an absorption line of 
fixed equivalent width varies as a function of S/N. The 'sig- 
nal' depends on the systematic change in the height of the 
quasar continuum+Clll] A1909 emission line with <5A. The 
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Figure 2. The results of testing for a significant difference be- 
tween the observed number of absorbers relative to the location 
of the Cm] A1909 emission line. The thick solid line shows \ 2 
for the pure data. The dotted solid line shows \ 2 once the ob- 
servations are corrected for the partial coverage effect (equation 
O; this results in more absorbers entering our sample, changing 
X 2 ■ The original data is consistent with predictions, unlike the 
corrected data, so that we conclude the correction is fictitious. 

'noise' varies due to a combination of the increasing sky- 
background as a function of observed-frame wavelength and 
of the sensitivity of the SDSS spectrograph and detector. By 
combining these considerations we can predict the number 
density of absorbers as a function of 5X. 

3.2 Statistical Test 

We adopt a very simple statistical test, comparing the pre- 
dicted number of absorbers with Wq 796 JS0.9A against 
the observed number, yielding a value of \ 2 ~ (M>in — 
N cxp (5\)) 2 /N C x P (5\) in bins of 40 A, where jV bln is the ob- 
served number of absorbers in the bin and N exp (SX) is the 
expected number of absorbers in the bin according to the 
procedure detailed above. The bin size is chosen to maxi- 
mize the number of absorbers per bin whilst ensuring the 
effect of the broad emission line is not diluted by choosing 
the bin width too large. The central bin is centred on S\ — 0. 

Since one expects, in the absence of partial coverage, 
iVbin to be Poisson distributed with mean jV exp , and iVbm ~ 
300 S> 1, the distribution should be accurately Gaussian and 
\ 2 should be ~ 1 in all bins. However, in the small absorber 
scenario, the dilution effect will cause a drop in the number 
of observed absorbers in the vicinity of the broad emission 
line, and hence an increase in x 2 ■ 

The result is shown by the thick, solid line in Figure [2] 
The mean \ 2 — 0.7 indicates a result consistent with the 
full coverage picture. To investigate whether this could oc- 
cur simply because the effects we look for are too small, we 
now perform a second test, in which we make the assump- 



tion that the observed equivalent widths are altered over 
the Cm] emission line due to partial coverage. In this pic- 
ture, we have observationally measured W' , so each absorber 
is then adjusted according to ([5| to yield W. C is deter- 
mined by interpolating continuum points from each side of 
the C ill] feature, whilst L+C and hence L may be estimated 
by inspecting the signal to either side of the Mgn feature. 
This results in a number of absorbers, which previously had 
equivalent widths below the 0.9 A limit, entering our sam- 
ple. The correction factor to equivalent widths of absorbers 
on top of the emission feature peaks at 1 + L/C ~ 1.5, al- 
though there are a few outliers where this factor takes a 
value as high as 1 + L/C ~ 2.1. 

Performing the statistical test yields an exceptionally 
poor fit between the observations and the model predic- 
tions, shown by the dashed line in Figure [2] We conclude 
that the partial coverage scenario is inconsistent with the 
observations at a very high level of significance, certainly 
in the simple limit considered here. We now show that the 
comparison x 2 results are, in fact, not invalidated by our 
oversimplification of the absorber model. 

3.3 Complication A: Off-Centre Single Absorber 

Suppose the continuum region has area A c , the absorber 
A a ^ Ac and the broad emission line region has area Ai. 
If the area overlap between the absorber and continuum is 
Ao, with total absorber area A a , it is simple to show that 
the observed equivalent width W' is given in terms of the 
intrinsic equivalent width W by: 

w > ApAjC + (Ag — Ap)A c L 

AoMC + L) ■ (<> 

Since one could never observe the intrinsic W in this pic- 
ture, the correct quantity to consider is W'(L, Aq) /W' {L = 
0,Aq>), i.e. the ratio between the equivalent width observed 
over an emission line to that observed where there is no 
emission from the line region: 

W'(L, A ) = ApAiC+jAa- A )A C L 
W'(L = 0,Ao) (C + L)A A l ' [ > 

For A = A a , so that W'(L = 0,^o) = W, this equa- 
tion correctly reduces to equation ([5]). However its behaviour 
elsewhere is somewhat subtle; in the limit Aq — > 0, one ob- 
tains W'(L, Ao)/W'(L = 0,Ao) — > co, apparently compen- 
sating for our original effect and generating extra high W 
absorbers. However, inserting the advocated A c > 4A a , and 
the observational values Ai > lO^c, L ~ 0.2C, one can use 
(0 to show that the population of absorbers generating this 
compensating effect must have intrinsic W > 250 W' . Thus 
those absorbers ever to enter our sample (W' > 0.5A) would 
have physically absurd high intrinsic column densities. We 
conclude the effects are only significant for W' <C 0.5 A, so 
that our statistical result (Section l3.2|) is insensitive to these 
considerations. 

3.4 Complication B: Multiple Absorbers 

It is clear from the multi-component nature of observed 
Mgn systems that any single cloud approach is over- 
simplified. If one asserts the overall size of the multi-cloud 
complex is smaller than the continuum region, our result 
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holds identically since we have simply shown that there is 
no detectable net difference in the fraction of light absorbed 
against the continuum and C ill] line regions. 

A separate scenario is multiple clouds per sightline, each 
with A a ~ A c . This would introduce absorption over the 
broad emission li ne, thus reduc i ng th e magnitude of our 
comparison \ 2 - In lPorciani et all ((2007) , it is stated that the 
number of clouds will scale with the beam area. This may 
or may not be a good approximation to the true behaviour; 
however, in our opinion the effect cannot be addressed in 
detail without a sophisticated model which includes con- 
siderations of not only the abundance but necessarily the 
spatial correlation of Mg II systems. 

For instance, in a model where Mgll systems cluster 
in spheres of a characteristic radius, the total number of 
absorbers over the line region will scale as r 2 /A c where r 
is the radius of the cluster of absorbers, for Al > trr 2 > 
A c , but as Al/A c for Tvr 2 > Al- This latter scaling is also 
applicable in the case that Mgll systems are considered to 
be evenly scattered across the local patch of sky. 

Discussion of these models are beyond the scope of this 
letter. But such considerations do not affect our conclusions 
with regard to the small Mgll cloud models. The mod- 
els of F06 are calculated considering only isolated clouds, 
and already these clouds must be adjusted in size to fit the 
data. In order to simultaneously fit our negative result and 
the GRB/quasar mismatch, one would need additionally to 
model the distribution of clouds so that the approximation 
breaks down in the correct manner and at the correct scale 
to compensate precisely for the geometric dilution in our 
specific case. Such a degree of fine tuning would not be 
acceptable without the support of a physically motivated 
model. 



4 CONCLUSIONS 

We investigated the possibility that Mgll absorber clouds 
are somewhat smaller than, but comparable to, the beam 
width of the quasar continuum. Such a scenario is motivated 
by disparate statistics of Mgll absorption when comparing 
GRB and quasar sightlines (F06). 

Ah our results suggest that, without extreme fine- 
tuning, the model is hard to match with the known proper- 
ties of Mgll systems. High resolution spectroscopy directly 
rules out partial coverage clouds with hard edges, whilst 
softening the edge will mimic a pure Voigt profile well but 
introduce difficulties with producing the known equivalent 
width distribution. Most convincingly, a statistical test for 
systematically lowered equivalent widths over quasar broad 
line emission regions (which are substantially larger than 
quasar continuum regions) produces no result. After apply- 
ing a correction on the assumption that partial coverage is 
the correct scenario, an inconsistency between the equiva- 
lent width distributions over broad emission line and con- 
tinuum regions arises. The only circumstances in which the 
F06 model holds and our broad emission line effect might 
be unobservable are those in which extreme fine tuning of 
the sizes and distributions are invoked (Section 13.41) . 

We therefore reject the proposition of Mgll absorbers 
which are patchy on scales comparable to quasar continuum 



beam sizes. The tru e explanation for the observations of 
iProchter et alj (|2006h remains elusive. 
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